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Abstract 

We evaluate analytically the elastic vrvr scattering amplitude to two loops 
in chiral perturbation theory and give numerical values for the two S'-wave 
scattering lengths and for the phase shift difference 6q — Sl- 



1. In the framework of chiral perturbation theory (CHPT) |1|, 0, the elastic vrvr 
scattering amplitude is evaluated by an expansion in powers of the external momenta 
and of the light quark masses, 

A = A2 + Ai + AQ + --- , (1) 

where An is of 0(p"). The first two terms in this expansion have been extensively 
analyzed during the last three decades 0, ^, Q . After a long period in which hardly 
any data have been collected, tttt scattering will receive in the near future interesting 
new input from the experimental side: i) It is expected that forthcoming precise data 
on Kei decay at DA$NE and at Brookhaven will allow one to determine the phase 
shift difference ~ "^i near threshold with considerably higher precision [Q than 
hitherto available ||10||. ii) There are plans to measure the lifetime of the vr+Tr^-atom 
in the ground state ||ll| , and to accurately determine in this manner the combination 
Oq — Og of the two S'-wave scattering lengths. In order to confront these data with 
precise theoretical predictions, it is necessary to go beyond the next-to-leading order 
term ^4 |^ . As has been pointed out in [^] and in [|I^ , one may eventually obtain 
experimental information on the size of the quark-antiquark condensate in QCD in 
this manner. 



In Ref. fT^, the part of the amplitude Aq containing branch points - required 



by unitarity - has been determined. A general crossing symmetric polynomial of 
O(p^), not fixed by unitarity, has been added by hand. In this letter, we present 
the analytical result for Aq based on a full two-loop calculation in the framework of 
CHPT. We compare the results of the two approaches below. 

2. The expansion (l|) is most conveniently performed in the framework of an effective 
lagrangian 0, ^|. Here we consider an expansion around the chiral limit m„ = 
irtd = 0, whereas the strange quark mass is kept at its physical value. We ignore 
isospin breaking effects and put m„ = = ih. The effective lagrangian is expressed 
in terms of the pion field U and of the quark mass matrix x, 

= C2{U, x) + hC,{U, x) + h'Ceiu, x) + • • ■ , 

where contains mi derivatives of the pion field and m2 powers of the quark mass 
matrix, with mi + 2m2 = n. Given , it is straightforward to expand the ^-matrix 
elements in powers of h. This procedure automatically generates the series (|l]), viz., 
A = A2 + hA^ + h^Ag + . . . . The leading-order term A2 has been evaluated by 
Weinberg whereas the next-to-leading order correction A4 was presented in 0. 
The calculation of Aq requires the evaluation of two-loop graphs with C2, one-loop 
graphs with one vertex from £4, and tree graphs generated by C2 + ^^.4 + Cq. Details 
of this calculation will be presented elsewhere. In particular, we refer the reader for 
the explicit expressions of £2,4 and of Cq to Ref. and Ref. , respectively. 



3. We use the notation 



(vr'^(p4)vr^(P3) out\TT''{pi)7T\p2) m) = {7t'^{pa)7t%P3) in|7r'^(pi)7r^(p2) in) 
+i{27ry6\Pf - Pi) [6''^6^'^A{s, t, u) + permutations} 
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where s, t, u are the usual Mandelstam variables, expressed in units of the physical 
pion mass squared M^, 

s ^ {pi+ P2?/Ml ,t={ps- Pif/Ml ,u={p,- pif/Ml . 

Using these dimensionless quantities, the momentum expansion of the amplitude 
amounts to a Taylor series in 



X2 = 



7r_ 

F3 



where F^^ denotes the physical pion decay constant. We find 
A{s,t,u) = X2[s-1] 



u 



with 
F(2)(s) 



+X2^ 



F^^\s)+G^^\s,t) + G^^\s,u) 
hs^ + b6s{t - uf] 
F^^\s)+G^^\s,t)+G^^\s,u) 
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+ ]-K3{t){l + 3st-s + 3t^-9t} + lK4{t) {4-2s-t} . (3) 
9 ^ ^3 

The loop functions J and Ki are displayed in appendix A, whereas the coefficients 
bi in the polynomial part are given in appendix B. 

4. We comment on the structure of the result. 

i) The amplitude A{s,t,u) is expressed in terms of the external momenta, the 
physical pion mass, the physical pion decay constant, and the coefficients 
bi,...,bQ. To arrive at this result, one has to evaluate also and to 
two loopsQ. Quantum field theory leads to the relations (|B.1|) , that determine 
bi in terms of 

1 

- chiral logarithms L = log — ^ , 

- the low-energy couplings /^(/i), . . . , IKfi) from £4 , 

- the low-energy couplings r[(yu), . . . , rg(/i) generated by . 

The scale dependent couplings /[ (r[) are needed to remove the ultraviolet 
divergences at order (p^), and to generate the most general solution of 
the Ward identities at these orders [0, The scale fi - introduced by the 
renormalization procedure - drops out in the full result (0). 

ii) The contributions proportional to Xg" in (□) correspond to terms of 0(p^"). 
The terms proportional to X2 in G^'^^ and generate contributions of 0{p^) 
- these are beyond the accuracy we aim at here. In order to keep the formulae 
as simple as possible, we nevertheless retain them. 

iii) We compare the result (|^) with the amplitude given in Ref. ||T3. Identifying 
the low-energy couplings a, /3, Ai, . . . , A4 introduced there with the relevant 
linear combinations of 61, . . . , fog) the two expressions agree at O(p^). The S- 
matrix method used in [jl3| and the field theory calculation presented here 



therefore agree as far as the absorptive part of the amplitude and the general 
structure of the real part is concerned. Our use of an off-shell method provides 
the additional information Eq. ( |B.1| ). Together with an estimate of the new 
couplings at O(p^), this allows us to make predictions e.g. for the S'-waves, see 
below. In Ref. [|13|, a and f3 are on the other hand treated as phenomenolog- 
ical parameters that appear in the expressions for the S'-waves which remain 
undetermined. Once a and f3 have been pinned down, that approach will even- 
tually allow one to determine the size of the quark-antiquark condensate, and 



to compare the result with the Gell-Mann-Oakes-Renner framework ||T6 | 



^ We are indebted to Urs Biirgi for providing us with the relevant expression for M^r prior to 
pubhcation ||l5||. 
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5. To proceed further, we need to know the values of the coefficients 6j. The 
constants /[ that occur in these have been determined from experimental data and 
using the Zweig rule some time ago [0] (for an update, see ^^). All these determina- 



tions are faced with the problem that the couplings /[ are quark mass independent, 
whereas the physical quantities, from where the are evaluated, incorporate quark 
mass effects. Here we have, for the first time, a means to pin down these quark mass 
effects at leading order in an algebraically precise manner. In order to achieve this, 
we need an estimate of the couplings r[. This is not an easy task, and we postpone 
a complete discussion to a later publication. Meanwhile, we use for /[ the values 
found in g 0, 

ll(Mp) = -5.40 • 10-3 , /^(Mp) = 5.67- lO^^ , 

lllMp)= 0.82-10-3, ll{Mp)= 5.60-10-3; = 770 MeV . 

As for r[, we use a by now standard method to obtain an order of magnitude 
estimate: We incorporate in the above representation the contributions from the 
lowest heavy states and assume that these effects account for the bulk part in the 
low-energy couplings at 0{p^). This procedure works very well at O(p^) P, |l^. In 
particular, we include the contribution from vector and scalar exchange, using the 



couplings presented in Ref. |T8[. In addition to these, kaons and etas also generate 
contributions of 0{p^). To estimate those, we have taken from Ref. |jl3| (see also 
Ref. |]19[) the elastic vrvr scattering amplitude of O(p^) evaluated in the framework 
of SIJ(3) X SU{3) and restricted it to SU{2) x SU{2) by an expansion in inverse 
powers of the strange quark mass. 



6. It is now straightforward to extract the scattering lengths af and slope parameters 
h\. For the isospin zero S-wave, we obtain 
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and similarly for the other threshold parameters. Numerically, we findQ by keeping 
terms up to and including 0(2:^3) 



0.217 (0.215) 
0.258 (0.256) 



(4) 



where the numbers in brackets denote the values obtained by putting the couplings 
at O(p^) to zero at the scale \i= 1 GeV, r\{\ GeV) =0. It is seen that for these 
threshold parameters the new couplings of O(p^) contribute a negligible amount, if 
their values are estimated in the manner described above. It is also worth empha- 
sizing that there is essentially no scale dependence for 0.5 GeV < < 1 GeV. We 



^For ease of comparison with earlier calculations, we use F-^ 
MeV. See also point 9. below. 



93.2 MeV 27fl and = 139.57 
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comment on the theoretical uncertainties of these and other predictions below. The 
result @) should be confronted with the data 



0.26 ±0.05 , 
0.29 ±0.04 . 



The value for is from Ref. 



2g 



21 



whereas we have used the universal curve 

to express Cq through and to obtain the second relation. Finally, we recall the 
result for the one-loop approximation 

= 0.201 (one-loop result) , 
O'o — clq = 0.242 (one-loop result) . 



7. The chiral expansion of the S-wave threshold parameters contains chiral loga- 
rithms ||2^, ^. At one loop, these are responsible for the bulk part of the correction 
to the tree-level result, if the scale in the logarithm is taken at 1 GeV 0. At two 
loops, the expansion contains also squares of chiral logs ||23|, e.g.. 



7X2 

327r 



with 



X2 



+ X2 

±0(a;/) } 
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-L ± analytic 

58, 96, ^, 11, 1697 L , . ' 

—ki + —k2 + 5fc3 ± — fc4 ± ^"TT^ + analytic 
7 7 2 84 167r^ 



(5) 



kii^i) = (4/[(/i) - 7,L)L ; 71 = 3 



72 = 3 



73 



'2 '^^ 



2 . 



(6) 



The coefficients of the fcj's in the threshold parameters, in particular in Eq. (|^), have 
been evaluated earlier in Ref. p3[ by means of renormalization group techniques 
M. Evaluating the expressions at the scale /x =1 GeV gives 



tree 



1 loop 



2 loops 



a° = 0.156 ±0.039 
L 



0.005 ±0.013 ±0.003 ±0.001 
anal. ki L anal. 



total 
0^ 



0.201 ±0.036 ±0.006 ±0.012 ±0.003 ±0.001 = 0.258. 



anal. 



anal. 



We conclude that the nonanalytic terms also dominate the two-loop corrections in 
this case. 



8. In Fig. 1, we display the phase shift difference 5q—5\ (in degrees) as a function of 
the center of mass energy i^^^^ = Mtta/s. The dotted (dash-dotted) line stands for 
the Born term (one-loop result), whereas the solid line shows the result at two-loop 
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Figure 1: The phase shift difference 6q—61 (in degrees) as a function of the center of 
mass energy of the two incoming pions. The dotted (dash-dotted) hne displays the 
tree (one-loop) result, whereas the solid line denotes the two-loop approximation, 
evaluated with r[(l GeV) = 0. The data are from Ref. pl| . 

accuracy, evaluated with r[(l GeV) = 0. It is seen that the two-loop corrections 
are reasonably small also considerably above threshold. Using the couplings r[ 
estimated in the above described manner increases the two-loop result at 380 MeV 
by 0.4 degrees. 



9. In summary, we have presented the analytic expression for the elastic tctt scat- 
tering amplitude to two-loop accuracy. In order to estimate numerically the size 
of the two-loop contributions, we have determined the new couplings that occur 
at this order by saturating them with the polynomial contributions to the ampli- 
tude generated by heavy states. In the case of the threshold parameters and 
Oq, we then find that i) the new couplings are numerically negligible, ii) the bulk 
part of the correction is due to the presence of chiral logarithms, and iii) there is no 
manifestation of a strong enhancement of the two-loop contributions in this case. 

A more reliable numerical exploitation of the representation (Q), in particular the 
evaluation of the remaining uncertainties in the predictions, requires additional work: 
i) As we mentioned above, the quark mass effects in the determination of /[ must 
be investigated, ii) In view of the accuracy aimed at in future experiments [TT| , 



isospin breaking effects cannot be neglected any further. For example, electroweak 
radiative corrections must be properly taken into account for extracting [^]. To 
illustrate, using = 92.4 MeV p6[ instead of = 93.2 MeV increases the 



values (H) for aj] (oq — Oq) by 0.005 (07006). If we use the neutral pion mass instead of 
M^+, the effect goes in the opposite direction. In addition, real and virtual photon 
emission in the scattering process should be investigated. We defer these and related 
issues to a future publication. 
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A Loop functions 

Let 
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log 
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s 



N = IQtt' 



Using matrix notation, the loop functions used in the text are given by 



and 
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The functions s~^J and s~^Ki are analytic in the complex s-plane (cut along the 
positive real axis for s > 4), and they vanish as |s| tends to infinity. Their real and 
imaginary parts are continuous at s = 4. The combination NKi{s) is denoted by 
Ki{s) in Ref. P . 



B The coefficients 61, . . . , 6^ 



The quantities bi in Eqs. @ and (|]) stand for 
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(B.l) 



1 M2 

with L = log — ^, and where the ki are defined in Eq. (g). We have denoted by 

167r2 /i^ '-' 

/[ (r[) the renormahzed, quark mass independent couphngs from £4 (Cq), with 

/X— = . The scale dependences of r^ are fixed by the requirement /i— = 0. 

dfj, 167r2 dfi 
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